Superfluid equation of state of dilute composite bosons 
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We present an exact theory of the BEC-BCS crossover in the BEC regime, which treats explicitely 
dimers as made of two fermions. We apply our framework, at zero temperature, to the calculation 
of the equation of state. We find that, when expanding the chemical potential in powers of the 
density n up to the Lee-Huang- Yang order, proportional to n 3 ^ 2 , the result is identical to the one of 
elementary bosons in terms of the dimer-dimer scattering length am, the composite nature of the 
dimers appearing only in the next order term proportional to n 2 . 
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The BEC-BCS crossover first considered by Leggett 
[l| , and the recent experimental realization of Bose Ein- 
stein Condensates (BEC) of molecules made of fermionic 
atoms 0, S 0, @ have motivated a number of theoret- 
ical works. Indeed, thanks to Feshbach resonances, it 
is experimentally possible, with two fermions of mass m 
( 6 Li or 40 K) in different hyperfine states (we denote them 
as 'spin' j and J,), with scattering length a, to realize 
weakly bound molecules, or dimers, with binding energy 
Eb = 1/ma 2 (we take K = 1 in the paper). In particu- 
lar one may obtain a dilute condensate of molecules. A 
crucial quantity controling the physics of the condensate 
in this regime is the dimer-dimer scattering length czm- 
This is however a highly non-trivial quantity to calculate, 
since one has to solve a four-body problem to find it. In 
the case of a broad resonance, one finds aw = 0.6 a by 
solving the Schrodinger equation Q or resumming the 
diagrammatic series (7(. The study of a Bose-Einstein 
condensate of composite bosons, where all the theory is 
formulated in terms of fermions only, was started a long 
time ago d, B wi th steady progress. Quite recently, Pieri 
and Strinati [10( derived the Gross-Pitaevskii equation 
from the Bogoliubov-de Gennes equations. However, be- 
cause of their approximate scheme, they end up with the 
Born approximation 2a for the dimer scattering length 
am instead of the exact result. 

In this paper, we present an exact fermionic theory of 
a BEC superfluid of composite bosons in the low density 
range. Our framework is completely general. Our present 
work is a first step toward going to higher orders, which 
will be clearly more complex to handle. Here we restrict 
ourselves to the T = thermodynamics. We obtain for 
the expansion of the chemical potential 11 of our fermions 
of single spin density n in the BEC regime: 
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Except for the obvious first term (which implies /i < 0), 
this is exactly the result found, for fiBose = 2/z, by 
Lee, Huang and Yang (LHY) [nj for elementary bosons 
with density n, mass tub = 2m and scattering length 
o-b = clm- The identity of the mean field term is natu- 



rally somewhat expected. However, even if it is reason- 
able to expect in our case a corrrection of the LHY type, 
it is not at all obvious that the coefficient is the same. 
We will see that our derivation is quite involved and has 
no systematic mapping on a purely bosonic formulation. 
In other terms one expects the composite nature of our 
bosons to enter at some stage in the theory. We find in- 
deed that this happens, but only at the level of the n 2 
term in Eq.([T]). Hence we prove that, for our composite 
bosons, the LHY term itself is unchanged with respect to 
elementary bosons 12 1. 

In order to perform a low density expansion, we need 
a "small parameter" in our theory. The most conve- 
nient one turns out to be the anomalous self-energy A(fc) 
which, together with the anomalous (or off-diagonal) 
Green's function F(k), is the hallmark of the superfluid 
state in the diagrammatic technique [3] • We will indeed 
see that at low density A(k) is of order n 1 / 2 , which could 
be anticipated from the standard BCS calculation [jj, LSI] • 
Hence by performing an expansion in powers of A(fc) in 
Feynman diagrams, we actually perform a low density ex- 
pansion. The full Green's function G(k) and self-energies 
are related by the completely general standard equations: 

G(k) = g (k) + F(k)A*(k)g (k) (2) 

F(k) = -G(k)A(k)g Q {-k) (3) 



where we have set [Q (k)] 



S(jfc), with E(fc) 
k 2 1 2m + n and 



the normal self-energy, G 1 (k) = uj 
k = {k,to}. 

We proceed in a natural way by finding the expansion 
of the Green's function G and F in powers of A at fixed 
fx. The single spin density gives the " number equation" : 
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At zeroth order in A(fe) 



the result is obviously n = since, without condensate, 
there are no fermions at T = 0, fi < 0. From particle con- 
servation, the lowest order is given by the second order 
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term: 
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T 3 (k,k;k) [Go(k)]' 



(5) 



where TV depicted in Figfjjb), has been discussed in 
Ref. 0, E3] an d contains all the normal state diagrams 
describing the scattering of a single atom by a dimcr 
(actually in the involved vacuum Green's functions we 
have to shift the frequencies by the chemical potential 
fi). This includes in particular a term —Go(—k) which 
is just the Born approximation for T3. In writing Eq.([5]) 
we have made use of the fact that, at this order, the k 
dependence of A(fc) can be neglected as will be shown 
below, and we have just denoted the resulting constant 
by A. The frequency integral in Eq.© can be calculated 
by closing the contour in the upper-half complex plane 
Imw > 0, where Ga(k) is analytic. It can be proved that, 
except for the Born term, T^(k,k;k) is also analytic in 
this half- plane. Hence the only contribution in Eq.© 
comes from the Born term. However this term is the 
only one considered in the standard BCS theory on this 
BEC side. We end up with the very surprising conclu- 
sion that, at this order, all the detailed physics involved 
in the atom-dimer scattering is irrelevant and that the 
result is merely given by the standard BCS calculation, 
namely n 2 = m 2 \A\ 2 (2m\fi\] 1 ^ 2 ). This shows that 
A is indeed of order n 1 / 2 . 

We consider now the anomalous self-energy A(fc), in 
order to obtain our equivalent of the " gap equation" [l2| . 
A(fc) describes two atoms (k \,—k J,) which go in the con- 
densate. Quite generally the contributions to A(fc) are di- 
vided in two classes, so we have A(fc) = Si (k)+§2(k). The 
first class, the only one found in BCS theory, gathers dia- 
grams where these two fermions first interact through the 
bare two-body potential, with Fourier transform V(q), 
the second class containing all the other possibilities. In 
full generality the contribution of the first class, shown 
diagrammatically in FigfIJa), can be written: 



*i(k) = -^Vck-kijFffci) 



(6) 



from the very definition of the full Green's function F. 
Note that Si (k) is independent of to and, since the poten- 
tial is very short-ranged, it depends on k only for very 
high momenta. 

In the second class, where the two incoming fermions 
(k |, — fc I) do not first interact, we proceed to a A ex- 
pansion. The first order term is already included in S\ 
and particle conservation implies that the next order con- 
tains A A* A, the rest of the diagrams being made only 
of any number of normal state propagators Go(fc) with 
any number of interactions, as shown schematically in 
FigHlc) ■ Moreover these diagrams can not contain loops 
of normal state propagators, since this would correspond, 
in time representation, to the creation of particle-hole 



pairs. Such processes are impossible in the normal state 
at T = and /i < 0, where the free particle propagator 
is retarded. When these constraints are taken into ac- 
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FIG. 1: (a) BCS-like contribution 8i(k) (b) Structure of the 
lowest order normal self-energy (c) The diagram for 82(h). 

count, including the 'no first interaction' condition, one 
ends up with the conclusion that these normal state dia- 
grams have exactly been considered in Ref. 0] (with again 
a trivial shift of all the frequencies by //, as for Go(fc)), 
and denoted by <&, except for a subtle point which we dis- 
cuss below and is accounted for by the slightly different 
notation Hence: 



<5 2 (fc) = i|A| 2 A<I>'(fc,-fc;0,0) 



(7) 



In writing Eq.0 we have taken advantage that, to lowest 
order (see Eq.([6])), A(A;) is a constant independent of 
k. Hence in this third order term, we can take A(fc) as 
constant. On the other hand it is clear from Eq.© itself 
that A(fc) depends in general on k. The factor 1/2 is 
required to avoid double counting which arises from the 
presence of two factors A. 

The difference between $ and $' stems from the fact 
that <& is reducible, while <&' is not since it is a contribu- 
tion to the anomalous self-energy. Specifically $ contains 
the contribution — Go(fc)Go(— k) (this is the Born term) 
and also a term arising from the normal self-energy E(fc). 
Hence, in order to obtain one has to subtract from 
<!> these reducible diagrams. However exactly these same 
reducible diagrams appear automatically if we write from 
Eq. (|2l3p the series expansion for G " 1 (fc)F(fc)GQ 1 (— k) in 
terms of the (irreducible) self-energies A(&) and E(fc). 
Hence it is more convenient to add these reducible con- 
tributions on both sides of the equation for A(fc), in 
which case we have — Go(k)~ 1 F(k)GQ 1 (— k) in the left- 
hand side and $ appears in the right-hand side, instead 
of <&' (note that this manipulation is actually valid to any 
order in our expansion). This leads to: 

-F(*) = Go(fc)*i(k)Go(-fc) 

+ \ I A| 2 A Go(fc)G (-fc)<f>(fc, -k; 0, 0) (8) 

We then eliminate F(k) in favor of S\ by making use of 
Eq.®. The summation of the last term over k intro- 
duces the dimer-dimer scattering vertex T4(0,0;0) = 



3 



J2k Go(k)Go(— fc)$(fc, — k; 0, 0) evaluated at zero dimer 
energy. It is directly related Q to the dimer scatter- 
ing length by (87r/aTO 2 ) 2 T 4 (0, 0; 0) = <±na M /m. The last 
step in our procedure is the standard elimination of the 
interaction potential in favor of the scattering amplitude 
[lit . I n our casc tms quantity has to be evaluated at 
the energy fi, because our shift in frequency. After this 
step, <$i(k) can be taken as a constant Si, since all the 
momentum integrals are rapidly convergent. We end up 
with: 
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where we have simplified by S\ and made 5i ~ A in the 
right-hand side. When we substitute for |A| 2 its lowest 
order expression found above in terms of n 2 , we find for /i 
the mean field part of Eq. (fTJ) , with the appropriate dimer 
scattering length clm- 

The above is only the first step in our derivation. The 
natural continuation would be to go to next order in A, 
i.e. to order A 4 in Eq.® and order A 5 in Eq.©. This 
would lead to a contribution of order n 2 in Eq. (JTJ) . How- 
ever this expansion is not regular, as it would be the case 
if we had a gap between the ground state and the first 
excited state. Indeed there is, in our neutral superfiuid, 
a branch of the excitation spectrum which goes to zero 
energy when momentum is zero. This is the collective 
mode, physically identical to sound waves in the low en- 
ergy range, which is known as the Bogoliubov mode for 
Bose-Einstein condensates of elementary bosons. Natu- 
rally its existence is a fundamental property of the con- 
densate [l6[ . In the following we include only the contri- 
butions coming from this low energy collective mode. 

The propagator of this collective mode is a two-particle 
vertex and it is the generalization to the superfiuid state 
of the normal state dimer propagator T2(P). It enters 
our formalism in the following way. In Fig.l b) and c), 
the terms A* and A act as 'source' and 'sink' of fermions. 
They are required since, at T = 0, /i < 0, no fermions are 
present except coming from the superfiuid. However we 
can in general well think of having a dimer propagator 
going from A to A* (and replacing them) as shown in 
Fig. 3. This plays the same role for 'source' and 'sink', and 
gives diagrams which must be considered. It is easy to see 
that, in the normal state, they give a zero contribution 
(since there are no dimers in the normal state). But 
in the superfiuid state, this dimer propagator has to be 
replaced by the collective mode and the result is non zero. 
The terms we have to retain are just the modifications, 
with respect to the previous results, coming from this 
substitution. Actually we do not proceed immediately to 
a A expansion and our procedure is equivalent to series 
resummation to avoid singularities. 

To proceed we have first to find in our framework the 
collective mode propagator, more specifically in the low 
energy, low momentum range. It has a normal part T(P) 



and an anomalous part T a (P), which both depend only 
on the total energy-momentum P = {P, fi}. We write for 
them the equivalent of Eq. (|2l3jl . i.e. the Bethe-Salpeter 
equations: 



r = t 2 

T a = T 2 _ 



t 2 Y irr r 



V\ r , — r a + To. 



r Q r 

irr 



(10) 

(ii) 



where we did not write explicitely the arguments which 
are P or — P: for instance T 2 stands for T 2 (P) and T 2 _ 
for T 2 (—P). The normal (T[ rr ) and anomalous ( Tf rr and 
Tf rr ) irreducible vertices are analogous to self-energies. 

Then we expand these irreducible vertices in powers 
of A. Again from particle conservation the lowest order 
terms are second order. The result for Eq. lfTtj]) is depicted 
in FigO The 'normal part' (i.e. without the A factors) of 
the irreducible vertices involves clearly the normal state 
dimer-dimer scattering vertex T 4 considered above, since 
all 'in' and 'out' lines are dimer lines. Again, at this 
lowest order, A can be taken as constant. In this way 
Eq. lflO]) and ([TTJ) become: 



T = T 2 



-T 2 |A| 2 T 4 r 

|A| 2 f 4 _r a 4 



T 2 A 2 T 4 T a (12) 



T 2 _ A* 2 T A T 



(13) 



where T 4 = T 4 (P/2, P/2; P/2), T A = (1/2)T 4 (P, 0; 0) and 
T 4 = (1/2)T 4 (0,P;0), the factor 1/2 being again topo- 
logical. 
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FIG. 2: Diagrammatic representation for Eq. (|12|l for T 

We can now solve for V and r a . In the low en- 
ergy limit |P| -C 1/a and <C 1/ma 2 , we find eas- 
ily r(P) = -8Tr(n + n B + P 2 /Am)/(m 2 aD) and T a (P) = 
87TfiB/(m 2 aD), where D = (P 2 /4m) 2 + 2^ B P 2 /4m-f} 2 . 
We have set /ib = |A| 2 maaM/4 and evaluated the factor 
of A to zeroth order by taking 2|/i| = 1/ma 2 . The collec- 
tive mode frequency is obtained by setting D — and we 
recover as expected the Bogoliubov dispersion relation. 

We now consider the additional contributions to the 
self-energies coming from the collective mode. For the 
normal self-energy we have to add the top diagram in 
Fig|3l which gives an additional contribution n cm to our 
lowest order result Eq.©: 
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T 3 (k,k;k + P)T(P)[G Q (k)} 2 (14) 



Actually we should have subtracted from T(P) its ze- 
roth and second order terms in the series expansion in 



4 



powers of A (this is indicated in Figj3] by the slash in 
the mode propagator), since they are in principle taken 
into account in Eq.©. However it is easily seen that 
they are zero since they contain normal state propa- 
gator loops. In Eq. (fT4")l we can first perform the in- 
tegration over the frequency variable u> of k, by clos- 
ing the contour in the upper half-plane. Just as above 
in Eq.Q, it can be proved that the only contribution 
comes from the Born term of T%{k, k;k + P). Then the 
k integration is easily performed and we are left with 
n cm = (m 3 / 2 /8^) Ep r(P)/[2| M | + P 2 /4m - fi] 1 / 2 . The 
n integration can be transformed over a contour which 
encloses all the singularities of T(P) on the real negative 
axis. The high energy contributions to n cm coming from 
1^1 ~ 1 /ma 2 (physically linked to breaking dimers) will 
give negligible regular terms of order A 4 , as discussed 
above. On the other hand the contribution of the low 
frequency collective mode is easily calculated with the 
low energy expression of T(P) given above. We find: 
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(15) 



where we have used the zeroth order expression 
for \fj,\. When we use for /j,b its lowest order expression 
in terms of n, we find that n cm coincide formally with 
the 'depletion of the condensate', known for elementary 
boson superfluids. 



to subtract from F a (P) the lower order terms already 
taken into account in our lowest order calculation, lead- 
ing to Eq.([9|). In contrast with the case of T(P), the 
term we subtract is not zero, but acts to regularize the 
remaining integral over momentum P, which would oth- 
erwise have a high momentum divergence 17|. We ob- 
tain finally for the slashed contribution, which takes into 
account this subtraction, f (P) = 3^ p r(P) = 
247m cm (2m|//|) 1 /2/ m 2 

If we gather all the contributions, we have for the sin- 
gle spin density n — ri2 + n cm while the " gap equation" 
Eq.© is changed into the simple form: 



2 2 



a M |A| 2 + 5naa M n cm (16) 



When |A| 2 is eliminated between the "gap" and the 
" number" equations, the consistently expanded result for 
fx is indeed found to be Eq.{T]). 

In conclusion we have shown how an exact purely 
fermionic framework can be used in the BEC regime of 
the BEC-BCS crossover, and we have specifically demon- 
strated that the Lee-Huang- Yang result for the chemical 
potential remains valid for the corresponding composite 
bosons. We are very grateful to M. Yu. Kagan for stim- 
ulating discussions at the beginning of this work. 
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FIG. 3: Collective mode contributions to the self-energies 

We proceed now in the same way for the collec- 
tive mode contributions to the anomalous self-energy. 
Corresponding to the diagram FigQJ), we have to 
add the two bottom diagrams in FigO Just as 
in Eq.([7]) we should take only irreducible diagrams 
into account. But handling this problem in the 
same way by adding the reducible contributions on 
both sides of the equation, we end up with Eq.® 
except, in the right-hand side, for the additional 
contribution A G (k)G (—k)$(k, -k; 0, 0)^ T(P) + 
(1/2) J2p r °( p )]- As in E q-0 the factor 1/2 is topolog- 
ical. Then we follow the same procedure as after Eq.©. 
Finally, as in the calculation of n cm , we retain only in 
the summation over P the low energy contribution, the 
other ones giving higher order terms. The summation 
J2p F(-P) has already been found in the above calcula- 
tion of n cm . The summation T a (P) is more involved 
since, as we mentionned earlier below Eq. (|14p . we have 
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